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Abstract. New Integral Field Spectroscopy of the central region of NGC 7331 reveals strong Ha emission in the 
well-known CO and HI ring of NGC 7331. The [Nil] /Ha ratio indicates that a large scale stellar formation process 
is taken place at the ring in agreement with previous hypothesis about the exhaustion of gas in the inner to the 
ring region. The dynamics of stars and gas are not coupled. There is a ring of peculiar velocities in the ionized 
gas velocity map. These peculiar velocities can be well interpreted by the presence of an axisymmetric inflow of 
40 km/s at the inner boundary of the large-scale gaseous ring. We infer an inwards total flux of 1.6 MQyr -1 . This 
value is typical of the accretion rates in hypothetical large nuclear black holes. Despite the large differences in 
the scales of the nucleus and the gas ring of NGC 7331, we suggest that this inwards flux is feeding the nucleus. 
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1. Introduction 

Peculiar motions are frequently detected in active galax- 
ies where the existence of several kinematically distinct 
gaseous systems, some of them suffering radial outward 
movements, has been attributed to the influence of the 
k*" ■ active nucleus (Garcia-Lorenzo et al. 1999, 2001, Arribas 
'O ! et al. 1997, Mediavilla & Arribas 1993). The outflow of 
■ gas close to the nucleus is relatively easy to detect be- 
5^ ! cause of the large velocities involved (of the order of 100 
km/s) and of the direct illumination from the strong ac- 
tive source. Inflow of gas seems to be much more difficult 
to observe as it is presumably caused by less outstanding 
mechanisms than the nuclear activity. However, it is very 
important to detect since inflow should provide the mat- 
ter that will eventually accrete into the nuclear black hole. 
To study the influence of the environment on the nucleus 
we will present new results about the kinematics in the 
central region of NGC 7331. Some peculiarities previously 
detected in the ionized gas velocity field of this galaxy 
(see below) make it specially interesting for this kind of 
studies. 

The large scale gas distribution in NGC 7331 is ring- 
like as first detected by Bosma (1978) in HI, Telesco et 
al. (1982) from NIR colors and Young and Scoville (1982) 
in CO. The ring is confirmed in other wavelengths and is 
taken today as a prototype gas ring. 
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Bower et al. (1993) used optical long-slit spectra to 
show that their models without a central black hole fit 
the observational data. However, there is an unresolved 
LINER nucleus (see, for instance, Cowan et al., 1994), 
so this galaxy could harbour a massive black hole of 5 x 
10 8 M Q . This is supported by the motion of ionized gas 
([Nil] + Ha emission line) reported by Afanasiev et al. 
(1989) in the 0.2-0.4 kpc (2.8-5.6 arcsec) zone and by the 
discovery of a nuclear X-ray source by Stockdale et al. 
(1998) using ROSAT. 

The existence of a bar has been suggested from non- 
circular motions by Marcelin et al. (1994) and von Linden 
et al. (1996). However the I and K band photometry in 
the central regions carried out by Prada et al. (1996) did 
not indicate any significant barred morphology. 

In the central kpc, Bottema (1999) found that the 
emission line gas (for R < 40 arcsec) seems to rotate slower 
than the stars. An explanation of the observations would 
consist in an inclined and warped gaseous plane. 

Using 2D spectroscopy, Mediavilla et al. (1997) con- 
cluded that the kinematics of the stars and most of the 
ionized gas are decoupled. These authors found that the 
kinematic axes of the ionized gas velocity map are dis- 
torted and rotated with respect to the stars, something 
that can be interpreted in terms of radial movements. 
However, the reduced field of view of the available ve- 
locity maps (about 7" x 7" ) makes difficult to verify this 
interpretation. 
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In this letter we are going to present new Integral Field 
spectroscopy covering the central 30" x 30" of NGC 7331, 
with the aim of understanding the kinematics of the region 
surrounding the nucleus. 

2. Observations and reduction 

NGC 7331 was observed with INTEGRAL (Arribas et 
al. 1998) in combination with the fiber spectrograph 
WYFFOS (Bingham et al. 1994) the 15 of July of 2001. 
The spectral resolution and coverage were 4. 8 A and 1445 A 
respectively. The central wavelength was 6200A. From the 
wavelength calibration we estimate that the uncertainty in 
velocity determinations is better than 10 km/s. We used 
the SB3 bundle of INTEGRAL. This bundle consists of a 
central rectangle of 30" x 34" and an outer ring of 90" 
in diameter. Four exposures of 600s each were taken from 
the central region of NGC 7331. These frames were com- 
bined to increase the S/N ratio and to reject the cosmic 
rays. Posteriorly, the spectra were extracted, calibrated 
in wavelength, and corrected from throughput. The re- 
duction steps were made using the INTEGRAL data re- 
duction package running in IRAF. More details about the 
reduction procedure can be found in Arribas et al. (1991). 

For several spectral features of interest we have deter- 
mined a grid of values at the locations of the fibers at the 
focal plane and have interpolated a map of the selected 
feature (see, e.g., Arribas et al. 1999). 

In Fig. ^ we show the resulting intensity and veloc- 
ity maps. The stellar continuum corresponds to the ~ 
5730-6800A, spectral range. The stellar velocity map has 
been obtained by cross correlating the spectra in the ~ 
5600A-6850A wavelength interval using the nuclear spec- 
trum as template. We also present intensity and veloc- 
ity maps obtained from the multi-Gaussian fit of the 
Ha+[NII]AA6548, 6583 blend. The nuclear spectrum is 
very similar to the ones obtained from the same region by 
Filippenko and Sargent (1985) or Mediavilla et al. (1997) 
and shows a strong Ha absorption (Mediavilla et al. give a 
description of the nuclear spectrum of NGC 7331 in com- 
parison with that of M31). To avoid the influence of the 
strong Ha absorption in the determination of the veloc- 
ity and intensity maps, we have considered in the fits an 
extra component representing this absorption feature (see 
Mediavilla et al. 1997, for details). 

3. Observational results 

A very noticeable feature in Fig. [I] is the ring structure 
observed in the Ha intensity map. The two-dimensional 
distribution of the Ha emission consists of a nucleus, a 
region with very low emission and then the Ha ring that 
begins rather abruptly at about 7 arcsec measured in the 
east direction (~ 500 pc equivalent to a galactocentric 
radius of about 2 kpc, considering i — 75°). Therefore 
we are observing the inner part of the well-known large- 
scale gaseous ring. In contrast, the [Nil] emission decays 
continuously from the nucleus towards the outer parts 



and is very faint at the ring location. In the ring the 
[NII]A6583/Ha ratio is about 0.3, typical of ionization by 
a HH-like continuum. This indicates that the ring is repre- 
sentative of a large scale stellar formation process. On the 
contrary, the [NII]A6583/Ha ratio in the inner to the ring 
region is very high, about 5. This is typical of LINERS. 
At the nucleus we found a value [NII]A6583/Ha ~ 3 also 
typical of LINERS and similar to the one (2.7) reported 
by Mediavilla et al. (1997). It is also remarkable that the 
transition between the LINER and the starburst regions 
is very abrupt and takes place at the inner boundary of 
the Ha ring. 

Another important feature in the maps of Fig. [l] is the 
presence of a symmetrical distortion in the velocity maps 
corresponding to the ionized gas, also coincident with the 
inner boundary of the starburst ring. 

Leaving aside the S-shaped distortion present in the 
innermost region (see Mediavilla et al. 1997), the stellar 
velocity map in NGC 7331 is quite regular. However, the 
dynamics observed in Ha and [Nil] is very much distorted 
when compared with the pure rotation velocity field of 
the stellar system. Notice that the central 7" x 7" of the 
Ha map (Fig. ^) are in very good agreement with the Ha 
velocity map presented in Mediavilla et al. (1997). 

4. The large scale inflow 

The new map of Ha velocities (Fig. |l|) gives a global per- 
spective of the central region of NGC 7331 and can be 
interpreted in terms of a ring of peculiar velocities super- 
imposed to a normal rotation pattern. In Fig. |^ we repre- 
sent the velocity as a function of galactic azimuthal angle, 
v(8), along the ring of peculiar velocities for Ha and the 
stars. As it can be seen the Ha curve presents a global 
blueshift of 60 km/s with respect to the stellar curve (the 
[Nil] curve, not represented here, suffers also a blueshift 
of 23 km/s). After removing the global blueshift, the Ha 
curve is shifted by about 25 degrees with respect to the 
stellar one (see Fig. ||). The [Nil] curve corrected from 
blueshift is intermediate between that of the stars and Ha 
showing a shift of 15 degrees with respect to the stellar 
rotation. The amplitudes of the Ha and [Nil] curves are 
117 and 95 km/s, respectively. Thus, the kinematic shifts 
with respect to the stellar rotation imply a radial velocity 
component of about 55 km/s for Ha and 26 km/s for [Nil]. 
With an alternative method, measuring the departures of 
the velocity field from the systemic velocity along the mi- 
nor kinematic axis, we have obtained a similar value, 45 
km/s for Ha. 

The west part of this kinematically peculiar ring is 
relatively redder and the east-part relatively bluer. As the 
dust lanes in optical images clearly indicate that the west 
part is the closest, we can conclude that we are observing 
a contraction motion (see Fig. |^). Thus, we have detected 
gas flowing from the Ha ring to the central inner region 
with a very noticeable axisymmetry. The velocity is maxi- 
mum at the inner boundary of the large-scale gaseous ring. 
It then decreases towards the center, probably as V oc (R 
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Fig. 2. along the peculiar velocities ring. Triangles 
correspond to the stellar component. Circles joined by a 
solid line correspond to Ha. Circles joined by a dot line 
correspond to Ha corrected from blueshift (see text). 



/o) _1 , where V is the inward component of the gas veloc- 
ity at the inner boundary of the ring, R the galactocentric 
radius and p the density, due to continuity requirements. 

Interpretations other than axisymmetric inflow 
cannot be in principle disregarded to explain the 
radial motion. An elliptical motion would produce 
a departure from circular motion, without a net 
inflow. Also, the observed motions could be at- 
tributable to a pair of tightly wound spiral arms, 
i.e. streaming motions along the spiral arms due 
to a density wave could produce a similar observed 
velocity pattern. These effects are not axisymmet- 
ric, but as long as axisymmetry is not perfect, 
they should be considered as alternative interpre- 
tations. However, it would be more difficult to 
account for the observed global blueshift in the 
framework of these alternative explanations. In 
any event, we favor axisymmetric inflow as the 
simplest explanation. 

Regarding the global blueshift, it is noticeable that 
the magnitudes of the vertical (responsible of the global 
blueshift) and radial inward (causing the kinematic axis 
shifts) motions are in good quantitative agreement in both 
cases, Ha and [Nil]. This suggests that we are seeing gas 
flowing from the edge of the Ha starburst ring in all al- 
lowed directions. Thus, the global blueshift is explained 
by the same process that gives rise to the radial inward 
motion. 

There is good spatial coincidence between the peculiar 
velocities and Ha starburst rings (particularly good at the 
east region where the edge of the Ha ring is more sharply 
defined). Anyway, there is not significant differences in 




Fig. 3. A scheme of the nearly axisymmetric gas flow 
from the Ha ring to the central inner region detected in 
NGC7331. West is to the right and North is up. The in- 
ternal radius of the plotted ring is about 2 kpc. 



the results about the global blueshift and the radial in- 
ward motion at the galaxy plane if the velocity curves are 
determined along the Ha ring. The intermediate position 
of the [Nil] kinematics between the stars and Ha can be 
explained if more than one kinematic component is con- 
tributing to the [Nil] emission. This is supported by the 
existence of two-peaked [Nil] lines (Mediavilla et al. 1997). 

The inward velocity is much greater than the sound 
speed and a pressure gradient can be ruled out as the 
mechanism causing the gas inflow. However, it is of the 
order of the Alfven velocity if we consider a magnetic field 
of 10~ 5 Gauss. That means that the inflow of gas from 
the ring to the center could be induced by magnetic field 
pressure forces if a strong but in principle acceptable mag- 
netic field were associated to the ring (see Battaner et al. 
1988). Alternatively, the inflow can be also attributed to 
the presence of strong winds originated in starbursts. The 
existence of a large scale process of stellar formation in 
NGC 7331 was suggested by Young & Scoville (1982) to 
explain the presence of the CO ring in terms of exhaustion 
of the gas to form stars in the central region. The discov- 
ery of large scale stellar formation in the ring fits naturally 
in this hypothesis that is also supported by studies of the 
Mg2 index (del Burgo et al. 1997). 

The gas inflow from the ring could eventually be ac- 
creted into the central black hole. To do an estimation 
of the accreted mass rate, we consider that the inwards 
flux is 27rRHV/9. The gas density (p) is taken as 0.5 xlO -2 
M Q pc -3 (Israel & Baas 1999), the equivalent gas disc 
height as H = 467 pc (Bottema 1999) and the galacto- 
centric radius as R = 2 kpc. Therefore we estimate the 
infalling flux as 1.6 M©yr _1 of the order of magnitude of 
the accretion mass rate in a luminous seyfert nuclear black 
hole (Blandford et al. 1990). 
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For a low luminosity AGN, such as that in M81, the re- 
quired accretion flux has been estimated to be of the order 
of 5 x 10~ 5 M©/yr (Ho et al. 1997), assuming an efficiency 
of conversion between matter and energy of 0.1. As NGC 
7331 has a nucleus of comparable luminosity (Stockdale 
et al. 1998) it seems reasonable that a part of the radial 
macroscopic infall observed here can feed the NGC7331 
AGN. 

Considering the large difference between the scales of 
the nuclear accretion disc and the macroscopic ring ob- 
served here, the interpretation of the inward flux as the 
feeding of the central black hole can only be considered as 
tentative. Actually, this identification should be assessed 
under the conditions of continuity and stationarity. Even 
the inward velocity, as a function of radius, peaks at the 
inner ring boundary, it is not vanishing at lower radii. 
This fact favours the interpretation that the observed flux 
eventually reach the nucleus. 

5. Conclusions 

Some conclusions are worthy to be summarized: 

We have detected an Ha ring coincident with the inner 
part of the well known CO and HI ring. The [Nil] /Ha ratio 
at the ring is typical of HH-like ionization and indicates 
that a large scale stellar formation process is taken place at 
the ring. This is in agreement with the hypothesis of that 
the gas in the inner to the ring region has been exhausted 
by massive stellar formation, that continues at the ring. 
The physical conditions in the inner to the ring region are 
typical of LINERS and the transition between starburst 
and LINER takes place very sharply at the inner boundary 
of the Ha ring. 

Ionized gas peculiar velocities are found near to the 
inner boundary of the CO-HI ring, compatible with an 
axisymmetric morphology superimposed to the otherwise 
regular map of the ionized gas. The comparison with 
the regular stellar velocity map makes the assumption of 
an axisymmetric contraction very plausible. The infalling 
matter is observed at the inner boundary of the large-scale 
gaseous ring and likely continues towards the inner region. 
Although several mechanisms are possible to explain the 
origin of the inflow, the existence of strong winds at the 
ring wall fits naturally with the presence of massive stellar 
formation in the ring. 

Even if the accreting flux could have very different 
causes we would be observing infalling velocities of the 
order of 50 km/s and infalling gaseous fluxes of the order 
of 1 M Q yr~ 1 . This is higher than the required accretion 
rate for a central black hole. Probably we are observing the 
infalling motion feeding the black hole having its source 
as far as 2 kpc away from the center. 

The inwards flux axisymmetry seems to rule out, in 
this galaxy, mechanisms based on bars. Bars have been 
proposed as sources for feeding AGN's and this have been 
observed to operate in the central regions of some active 
galaxies (See Perez et al. 2000, and references therein). 
NGC 7331, without a clearly observed bar, and with an 



axisymmetric inflow suggest that other ways of loosing 
angular momentum could also be important. 
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